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MALYTICMi INVESTIGATION OF THE MOTION OF A 


CLAMSHELL-TYPE HEAT SHIELD ON DEPLOYMENT FROM A 

PARENT VEHICLE IN THE ATMOSPHERE 

By Ross L. Goble 
Langley Research Center 


SUMMARY 


The equations of motion for the deployment of the clamshell- type heat 
shield are developed herein for the various phases of separation from an accel- 
erating vehicle in the sensible atmosphere. The shield motion is planar and is 
divided into four phases: Phase 1, translation from the closed position to 

pivot contact; phase 2, rotation about pivot; phase 5^ translation on and rota- 
tion about pivot to release; and phase 4, translation and rotation away from 
the parent vehicle. Closed-form solution of the linear equation of motion in 
phase 1 for constant lift and drag permits description of the translation 
velocity. An iterative numerical integration technique is used to solve the 
nonlinear equations of angular motion for phases 2^ 5^ release equa- 

tion is also established. Specific applications of the analysis to deployment 
of the Nimbus heat-shield configurations^ used on many NASA projects, under 
various vehicular acceleration and aerodynamic loading conditions are included. 


INTRODUCTION 


Heat shielding is required for the protection of many types of flight pay- 
loads which have telemetry equipment, delicate instrumentation, and other 
reentry and space experiment apparatus. One specific requirement in the Project 
Fire experiment (ref. l), for example, was the shielding of the quartz windows 
in the forward section of the payload from atmospheric abrasion during the 
early flight stages. In some experiments, the shield is an integral component 
of the payload. Often, however, the nature of the experiment requires deploy- 
ment of the heat shield during flight to permit proper functioning of the pro- 
tected apparatus. Motion analysis of the shield intended for such use is 
required so that optimum design parameters and ejection conditions may be 
selected to ensure the desired deployment. The analysis herein is based on the 
spring method of shield deployment such as has been successfully used on the 
Nimbus shroud configuration, which includes the basic configuration used or 
intended for use on such projects as Fire, Echo II, Gemini Docking, Earth 
Orbiting Geophysical Observatory (eOGO), Nimbus, and other space experiments. 
Analysis also permits description of the shield trajectory after release. This 



description is important in determining the most desirable release conditions 
which will not only preclude shield mechanical interference with any portion 
of the vehicle after release^ but will also satisfy a predetermined minimum 
safe clearance requirement. 

The studies presented herein may be effectively utilized in the estab- 
lishing of optimim ejection criteria for future projects using the spring- 
ejected shroud. It is readily seen that many of the problems pertinent to the 
motion of the clamshell-type shield herein described will pertain to any 
clams he 11- shield separation regardless of the deployment means. 

It is convenient to consider the shield planar motion in four phases: 
phase 1, translation from the closed position to pivot contact; phase 2, rota- 
tion about pivot; phase 5^ translation on and rotation about pivot to release; 
and phase translation and rotation away from the parent vehicle. There are 
four parallel springs used in the Nimbus shield-ejection scheme, two acting at 
the aft end of the shield and two acting in the conical section near the nose. 
The springs act through phase 1 in translating the shield from the closed posi- 
tion to pivot contact. The two aft springs then serve the purpose of holding 
the shield pivot fitting against the booster pivot until shield release is 
effected at the end of phase 3. Also, the forward springs provide the force 
necessary to initiate and sustain rotation through phases 2 and The springs 
in this ejection scheme are designed to act until expanded to free length. 

Separation dynamics for vacuum flight have been studied for the two-body 
system in references 2 and 5 heat shields in references k- and 5- The 

consideration of aerodynamic effects has been generally confined to airplane 
escape capsules (ref. 6) and to atmospheric abort of space capsules from 
boosters (ref. 7)- The presented method is similar to that described in ref- 
erences 4 and 5 with the inclusion of aerodynamics. 

The purpose of this study is to investigate some of the factors which 
affect deployment and to ascertain their significance. The dynamic pressure, 
which has been generally neglected in related published literature, is one of 
the most important factors influencing shield deployment. Other factors treated 
in the application of the analytical method to physical systems are vehicular 
longitudinal acceleration, spring force, and change in length of the cylindrical 
portion of the shield. 


SYMBOLS 


Measurements for this investigation were taken in the U.S. Customary Sys- 
tem of Units. Equivalent values are indicated herein in the International 
System (Sl) in the interest of promoting use of this system in future NASA 
reports. Details concerning the use of SI, together with physical constants 
and conversion factors, are given in reference 8. 

A forward spring initial preload, in. (m) 

ai net vehicular longitudinal acceleration, g units 
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B 


aft spring initial preload, in. (m) 


Ci,C2,C3 

c 

D 

d 

Ek 

e 


®r 

Fa 

Fb 

Ffb 


H 

h 


^cg 


lO 


Kl 

k 

Fa 

Fb 

L 

L 


constants of Integration 

radial distance to shield center of gravity from pivot point, in. (m) 

aerodynamic drag force, Ihf (N) 

total shield diameter (fig. l), in. (m) 

kinetic energy, in-lbf (j) 

displacement of pivot along release plane, in. (m) 
displacement of pivot at time of release, in. (m) 
forward spring force (eq. (26)), Ihf (n) 
aft spring force (eq. (26)), Ibf (U) 
kinetic friction force at shield base, Ibf (n) 
static friction force at pivot, Ibf (n) 

acceleration equivalent to earth gravity, ^86 in/sec^ ( 9*8 m/sec^) 

pivot reaction in x-direction (fig. 2(b)), Ibf (N) 

distance from shield center of gravity to shield outside edge 
(fig. 1), in. (m) 

Integers 

half-shield mass moment of inertia about its center of gravity, 
Ib-sec^-in. (N-sec2-m) 

half-shield mass moment of inertia about the pivot center of rotation, 
lb-sec2-ln. (N-sec2_m) 

constant for phase 1 motion (eq. (8)), in/sec^ (m/sec^) 

superscript indicating kth Iteration for nimerical method, dimen- 
sionless 

forward spring constant, Ib/in. (N/m) 
aft spring constant, Ib/in. (N/m) 
aerodynamic lift force, Ibf (n) 

Lagrangian, Ek - Up,' in-lbf (j) 
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M 

m 

n,n+l 

P 

<1 


ll 




t 

Up 

V 

w 


Wp 

Wl 

% 

^,fb 


X 


Xp 


aerodynamic moment^ in-lbf (N-m) 

mass of shield half^ lbm-sec2/in. (kg) 

successive steps in the numerical solution, unitless 

distance from pivot center of rotation to contact point of aft 
spring, in. (m) 

dynamic pressure, Ibf/ft^ (w/m^) 

generalized coordinate having dimensions of coordinate used 

constant lased in clearance equations (k8) and (1^-9) 

distance from pivot center of rotation to contact point of forward 
spring, in. (m) 

distance from shield-half center of gravity to clearance point con- 
sidered (fig. in. (m) 

time variable, sec 

time intervals for deployment (fig. 9)^ sec 
potential energy, in-lbf (j) 

pivot reaction in y-direction (fig. 2(b)), Ibf (n) 

work done by a force, in-lbf (j) 

work due to aerodynamic drag, in-lbf (j) 

work due to aerodynamic lift, in-lbf (j) 

work due to aerodynamic moment, in-lbf (j) 

work due to friction force, in-lbf (j) 

longitudinal coordinate of shield center of gravity from reference 
axis (fig. 1), in. (m) 

longitudinal coordinate of clearance point from reference axis 
(eq. (49)), in. (m) 

pivot point longitudinal displacement used in virtual work (eq. ( 26 )), 
in. (m) 


location of shield center of gravity measured from shield base 
(fig. l), in. (m) 
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y lateral coordinate of shield center of gravity from reference axis 

(fig- 1), in. (m) 

y* location of shield center of gravity measured from shield parting 

plane edge (fig. 1), in. (m) 

Ya = y - ya> in- (m) 

ya = ®ro - cos(e + e^oj] ^ in- (m) 

ye = y - yh^ ^n. (m) 

yb = p(cos 9po - cos(e + 6po3> in- (m) 

y^^ shield center of gravity lateral displacement at pivot contact, 

in. (m) 


ye = yc - y* 



yp 


lateral coordinate of clearance point from reference axis (eq. (48)), 
in. (m) 

pivot point lateral displacement used in virtvial work (eq. ( 26 )), 
in. (m) 


yr = y “ y* 

A time interval prefix used in n\mierical solution 

6 prefix used with linear variables to indicate virtual displacement 

or virtiaal work 




0 

0o 

^co 

®po 

®ro 

M'S 


error terms used in numerical solution 

angle of rotation, shield base from original position, radian 
initial value of 0 for a given phase, radian 

reference angle between shield base and center-of-gravity location, 
radian (See fig. 2(a).) 

reference angle between shield base and aft spring point of con- 
tact, radian (See fig. 2(a).) 

reference angle between shield base and forward spring point of 
contact, radian (See fig. 2(a).) 

coefficient of static friction 

coefficient of sliding friction 
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I orientation of pivot relative to shield outer surface (fig. 2(b)), 

radian 

0j reference angle for jth clearance point (fig. 5)^ radian 

t = 0 - I, radian, (fig. 2(b)) 

Dots over variables Indicate differentiation with respect to time. 


DERIVAIEION OF EQUATIONS OF MOTION 


The derivation of the equations of motion for the rigid body is developed 
in this section for the four basic phases of operation. Also derived is the 
equation defining the conditions governing full release of the shield from the 
parent vehicle. Lagrange's equation is utilized in developing the equations 
of motion for all phases. 


Assumptions Pertinent to the Analysis 

The considerations -vdiich have been employed in the derivation in providing 
an adequate mathematical model to the physical system under investigation are 
as follows; 

(1) The vehicle is the inertial reference for the heat shield motion. 

( 2 ) At the inception of motion the starting friction force is neglected. 

( 3 ) The coefficient of kinetic friction is assumed to be independent of 
the translational velocity on the pivot. 

(4) The Coriolis effect dxiring translation off the pivot is neglected. 

( 5 ) The change in the center of rotation during translation off the pivot 
is considered to be sufficiently small as to be negligible. 

(6) The spin rate of the parent vehicle is assumed to be negligible 

in the analysis. It is also ass'umed that no roll motion of the shroud halves 
occurs during deployment. 

( 7 ) The heat shield is assumed to be a rigid body. 

(8) The two halves of the clamshell shield act independently. 

( 9 ) The aerodynamics including flow-field interference effects of the heat 
shield half are known. Induced aerodynamic loading due to deployment veloci- 
ties 0, y, and x have not been considered in the derivation. 
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Phase 1 - Translation to Pivot Contact 


If the coordinate system in figure 1 is used, the potential energy of the 
system is; 

Up = I ka(A - Jrf' + | " Vrf' (l) 

where y^. = y - y*, the displacement of the shield parting plane edge, and A 
and B are the preload deflections of the forward and aft springs , 
respectively. 

The kinetic energy is 



( 2 ) 


The Lagrangian L is given by - Up, 
system is 


and Lagrange ' s eq.Toat ion for the 


dt ^ Sy 5y 


(3) 


The forces (fig. l) acting on the shield are the friction force lift, 
drag, aerodynamic pitching moment, and the longitudinal acceleration force ajm. 
The virtml work can be described as 

8W = + &Wl (4) 


Since 


Ffb 



+ a^m 


) 


then 

(5) 

also 

5Wl = 31^ Sy (6) 

The aerodynamic moment and 
the forces in the x-direction 
do no work since the plane 
motion is confined to pure trans- 
lation in the y-direction hy 



Figxire 1.- Heat-shield model. Phase 1, translation. 
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constraints which preclude rotation and x-translation in this phase. Using the 
general relation for virtual work 


OO 



and expanding in terms of the single generalized coordinate y yields 


m = ~ by 
dy 

or^ after substitution of equations ( 5 ) and ( 6 ) for 6 W, 


-|ig(D + a^m)6y + L 6y = 8y 


Simplifying and equating coefficients^ therefore^ yields 

0 = L - Hs(d + a^m) 

Operating on equation (3) with substitution of the Lagrangian and 
yields 


y - kg , (a - y + y*) - k^(B - y + y*) = L - Hs(d + a^m) 


my 


Rearranging yields the differential equation 


y ( — in" ^) y = sl^aU + y*) + kb(B + y"") + L - Us(d + a^m)| 


Bw/Sy 

( 7 ) 


(8) 


Since motion in this phase is nearly impiolsive, it is assumed that aero- 
dynamic interference effects between the opening clamshell halves can be neg- 
lected. Thus, the aerodynamic forces in the absence of changing angle of attack 
may be considered to be constant for this phase. The right-hand side of equa- 
tion (8) is therefore a constant which may be called The complete solu- 

tion of the nonhomogeneous equation (8) is thus: 

y = Cl cos ^ - — ^ t + C2 sin ^ t C3 ( 9 ) 
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Substituting the particular integral y = Cj Into equation ( 8 ), with the 

right-hand side of equation ( 8 ) heing replaced hy Kp, and differentiating 
yields 

Km 


Substitution of the appropriate initial conditions^ t = 0 ^ y = y*^ into solu- 
tion (9) yields 


Cl = - 


+ Y 

^a ^ ^b 


Differentiating equation (9) with respect to tijne and substituting the initial 
conditions t = 0 ^ y = 0 yields C2 = 0 ; therefore, the complete solution 
with arbitrary constants evaluated and some rearrangements made is 


y = y - 


\ ,ra + 

cos\/ t + 


^a 


m 




(10) 


from which the solution for time t is 


y - 


Kpm 


i 


m 


+ ^b 


cos 


^a + ^b 


- 


Kim 


^a + ^bi 


( 11 ) 


The translation velocity relation obtained by differentiating equation (lO) is 


y 





or by substituting for t 


from equation (ll) and simplifying is 



(12) 
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Knowledge of the cent er-of -gravity displacement (for exang>le, y = Yc) which 
pivot contact occ\ars permits solution of equation (l2). When y^ is deter- 
mined, the initial value of the angular velocity for phase 2 can be established 
in the following manner. 


Phase 2 - Rotation About Pivot 

If it is assumed that no energy is lost in impact, all translational energy 
is converted to rotational energy about the pivot at the instant of pivot con- 
t act , or 



(13) 


Making substitutions for y^. 


and simplifying yields 



(14) 


The top and bottom springs have displaced an amount Ye == yc “ y* 
instant of pivot contact. Thus, the potential energy of the system taking the 
datum as shown in figure 2 is 


Up = a2mc sin(e + Geo) + | k^(A - Ya - Ye)^ + | - Y^ - Yef (l5) 


where 


ya = - cos(0 + Gj. 0^1 

yb = pjeos Gpo - cos(e + 0pQ)j| 


(16) 


and equation (l5) is valid for ya + ye = -^ yb ye - kinetic 

energy of the system is 


Ek = I 


(17) 
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(a) Geometry and release plane description, (l>) Detail of pivot point 

and release plane. 


Figure 2.- Heat-shield model. Phase 2, rotation. 

If the substitution for and y-^^ is made in equation {± 3 ) , the Lagrangian 

then becomes 

L = i Io9^ - a^mc sln(0 + 0co) - J kajl - r[cos - cos(0 + 0i.o)] - 


- I ^b|® - P|f°s 0po - cos(0 + 0poj] - y^‘ 


( 18 ) 


Operating on L with respect to the generalized coordinate 0 by Lagrange’s 
equation 

A.l^\ ^ M 

<it\ss/ " he he 

yields 

Iq© + ajmc cos(0 + ©co) - kajl - r [cos ©^o - cos(© + ©roQ “ y^^ sln(0 + ©ro) 
- k-5^ - p[cos ©po - cos(© + ©po)]] “ y^P sln(© + ©po) = (l9 


11 


For a virtxial angular displacement 60 the virtual work contributions of the 
aerodynamic forces (where D, and M are with respect to the shield center 

of gravity) are 


and 


6Wl = Lc|cos^0 + 9co) " cos^S + Gqo ■*" 

(20a) 

5Wj) = -Dcjsin^G + Qqq + 80 ) - sin 0 + 0co^ 

(20b) 

8% = M 80 


Expanding the double angles and making the small angle 
80 = 1, sin 80 = 80 yields, after simplification. 

approximations 

8Wl = Lc sin^0 + 0co^80 

(21a) 

8Wp = -Dc cos^0 + 0co]59 

(21b) 

also 


8% = M 80 

( 22 ) 


Again, employing the relation 


8 ¥ 


00 

V 3“ 

' I SJT ^“11 

i=l 


and substituting for the virtual work terms with respect to the generalized 
coordinate 9 yields 

|lc sin ^9 + 9co) “ ^ cos ^9 + 9 qo) + 50 ” ^ 

Equating coefficients and making substitution for Sw/^9 into equation (l9) 
yields, for the moment relative to the pivot point. 
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where lift, drag, and pitching moment vary with angle of attack 0 and include 
interference effects. Also, induced aerodynamic loading due to S, y, and x 
has been assumed to be trivial. Solution of equation (25) for 0 yields 



which is valid for 


& 




r cos 


) c 


P cos 


0*jPQ ■“ cos ^0 "J" 

0po - cos^0 + 


9ro)] 

0po)J 


- ^0 

- S 0 


(2h) 


Phase 2 initial conditions are as follows; At t = to^ 0=0, 0 = 0 q, 

L = Lo, D = Dof M = Mo, and the initial acceleration from equation (2il-) 
becomes 


""loir ~ 

+ LqC sin 0CO + Mo + k-^^B - ye)p sin 0p<^ (25) 

Solution of the angular displacement 0 and its time derivatives may be 
accomplished by stepwise numerical integration of equation (24) as shown in 
appendix A. 
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Developed in the next section are the relations required to analyze trans- 
lational motion on the pivot. 


Phase 5 - Release Action Equations 

The externally applied loads (fig. 2) are L, M, D, and a^mj the reac- 
tions are H and V; and the spring forces are Eg, and F-jj. The virtual work 
of the system is then 

6 W = -(d + ajm)5x + L 6y + M 60 + H 6xp + V 6yp + Fa + F^ dyg (26) 

where 6xp and 6yp are the respective x and y virtual displacements of 
the pivot point, and 6yj^ and 6y-g are the y virtual displacements of the 
forward and aft springs. 

The system kinetic energy related to the center of gravity of the shield 
half is 


Ek 




Icg9‘ 


(27) 


Since the system is originally restrained at the pivot, the following geometric 
relations in terms of the generalized coordinates x, y, and 9 may be written 

y = x*sin 9 + h(l - cos 0) + y<j (28) 


yp = y + h cos 0 - x*sin 0 

y;^ = y + c cos^0 + 9co) - ^ cos^0 

yp = y + c cos ^9 + - p cos^0 


®ro) > 
+ 0po) 


X = x*cos 9 + h sin 0 
Xp = X - h sin 0 - x*cos 9 


(29) 


(30) 

(51) 


Note that relations for y-p and Xp vanish identically hut are in a form 
necessary for establishing the virtual displacement of points actually con- 
strained against possible motion. This useful artifice is commonly employed in 
virtual work relations in order to permit definition of reactive forces at a 
point. The time derivatives of equations (28) and (50) are: 



I 


y = 0(x*cos 9 + h sin 0) 


y = o(x*cos 0 + h sin 0) + 0 (h cos 9 - x*sin 0) 


X = 9(h cos 9 - x*sin 9) 


it = 0(h cos 0 - x*sin 0) - 0 (x*cos 9 + h sin 0) (; 

Making substitutions into Lagrange's equation and differentiating yields 

my = L + V + Fa + F-b (; 

mx = -(d + a^m) + H (■ 

and 0 is redefined in terms of the V- and H-components by 

leg*©' = M + H(x*sin 0 - h cos 0) - V(h sin 0 + x*cos 0) + Faj^r sin^0 + 0 j.q^ 
- c sin (9 + ©co) + Fbjp sin ^9 + 0po^ - c sin ^9 + ©co^ 


If y and x from equations (33) and (35) are substituted into equations ( 36 ) 
and (37) 3.nd the terms rearranged^ the pivot reactions^ with Fa and Ft) 
expressed in terms of the respective spring preloads and displacements, are 

V = m|^9(x*cos 9 + h sin 9) + 9^(h cos 9 - x*sin 9)^ - L 

- ka^A - ye - rjeos ©j^q - cos (0 + ©ro) j 

- kb|B - yg - p|cos 0po - cos^0 + 9po)^ (59) 


H = mj^9(h cos 9 - x*sin 9) - 9^(x*cos 9 + h sin 9^ + D + a^m (^) 

The instantaneous orientation of the release plane is shown to he at some 
inclination angle \[r relative to the pivot reaction components, (See fig, 2.) 
The release plane will be partially determined by the physical configuration of 
the system. At release, the static equilibrium equation to be satisfied along 
this plane is 
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V sin iff - Ffp + H cos ijr = 0 


(in) 


where the static friction force is defined as 

Ffp = |i(V cos i|f - H sin \|r) 

Positive H and V have been chosen, as shown in figure 2(h). The equilihriim 
condition is, after substitution for Ffp in ecjuation (^l) and simplifying, 

V(sin - |i cos )|f) + H(cos ijf + p sin i|r) = 0 (^2) 


After substitutions for V and H from equations (59) (^O) are made, the 

release relation (42) becomes 



Equation (45) can thus be evaluated for any 0 and is the condition which 
establishes impending translational motion on the pivot. It can be shown that 
with slight modification, equation (45) becomes the condition of dynamic equi- 
librium on the pivot for the translational phase. This modification involves 
replacing the right-hand side of equation (45) with -me and n with |ig. 
Since the magnitude of the Coriolis effect due to the shield's rotating about 
while translating along the pivot is very small relative to other forces in the 
physical system, the slight retarding effect to the angtilar motion has been 
neglected. Also, for the small pivot translational displacement involved, the 
resultant change in the center of rotation can be neglected, and thus the geo- 
metric relations previously established as adequately describing the system are 
preserved. Therefore, the equation of translational motion (eq. (45) as modi- 
fied) is, after rearranging. 


l6 



^~|a - ye “ rj^cos - cos(e + 6^0)]^ 
- cos^e + 9po) I - |^(x* 


m 


B - 7e - p cos 9po 


'*’cos 9 + h sin 9 ) + 9^(h cos 9 - x*sin 9 




m , 


(sin \|f - jxg cos i|r) 


-([• 


9(h COS 0 - x*sin 0 ) - 9^( x*cos 0 


+ h sin 0) + D + a^m^(cos ^+ 1-13 sin i|r) 


m 


Solution of the translational displacement e and its time derivatives is 
accomplished by the numerical integration technique indicated in appendix B. 
When the shroud has translated on the pivot a distance e^^ release is effected 
and the fourth phase of motion begins. 


Phase k - Free Trajectory Subsequent to Release 

At release , the pivot reactions vanish, and the shield begins to translate 
and rotate away from the vehicle. The initial translational displacements and 
velocities for this phase are solved by equations ( 28 ), ( 30 ), ( 32 ), and ( 3 ^)^ 
respectively, evaluated at 0 eqtial to the release angle as determined at 
e = e^. Initial conditions for 0 q and 0 q are also those values established 
at release. The describing equations of motion, equations ( 36 ) to ( 38 ), become 

my = L (if5) 

mx = -(d + (46) 

and 

Icg9 = M (4?) 

Initial conditions for the accelerations are now readily obtained by substitu- 
tion of the proper values of the aero(3ynamic effects L, M, and D as deter- 
mined at 0 equal to the release angle. Solution of 0 and the corresponding 
y- and x-coordinates is accomplished numerically as outlined in appendix C. 

During the phase 4 motion as described, knowledge of the minimum clearance 
of the shroud halves from the booster vehicle is important, since any contact 
could cause damage to the booster and result in possible mission failure. By 
utilizing the geometry of the shroud, the minimum clearance envelope can be 
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described by relating coordinates of 
various points on the shroud to the 
coordinates of the shroud center of 
gravity as determined by solution of 
the displacements from eq.uations (^5) 
and (46). (See appendix C.) Fig- 
ure 3 indicates the six shroud loca- 
tions chosen to describe the clear- 
ance envelope. The geometric 
relations utilized are the following: 

yj = y + qjrj cos(e + (48) 

^ (^9) 

where 

qj = 1 (j = 1 an<i 2) 

qj = -1 (J = 3, 4, 5, and 6) 


The absolute displacement of the 

Figure 3.- Clearance point identification. shroud minimum point can readily be 

detennined by relation (48) at the 
end of each time interval. Rela- 
tion ( 49 ) can now be solved for the corresponding longitudinal displacement. 
The lateral clearance is then determined by differencing this near point 
yj -displacement with the booster outer radius at the x-coordinate . 


GENERAL DISCUSSION AND RESULTS 


The equations of motion for a clamshell spring-actuated heat shield^ 
deployed under conditions of aerodynamic and constant vehicular longitudinal 
acceleration loading, have been derived and presented in the foregoing sections. 

In the following sections are presented various factors involved in the 
application of the analytical equations to the Nimbus long and short heat- 
shield configurations. Results of a parametric study to determime separation 
behavior with variations in dynamic pressure, vehicilLar acceleration, spring 
force, and shield length are provided. 
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Simplified Shield Geometry 

Mecha nical considerations . - The clamshell- type heat shield as used on 
Project I^e~Is of the type for which the foregoing analysis is applicable. 
Several comments pertinent to the actual shield configiiration which verify cer- 
tain aspects of the analytical approach are as follows: First, during the 
phase 1 motion, the constraint which precludes rotation is provided by a nose 
latch mechanism. Although this constraint in reality permits a slight nega- 
tive rotation during this phase, the trajislational distance at the shield base 
is so small that the assumption of translation only is valid. Secondly, the 
spring act\iator system is such that for symmetrical shield halves , the springs 
operate in a direction always perpendicular to the vehicle longitudinal axis. 

Symmetrical halves .- As previously discussed, the assumption of symmetry 
in the shield-half geometry and spring force distribution was made to simplify 
the analysis. However, due to near symmetry, it was fo\md that either of the 
actual shield-half physical parameters such as inertial and geometrical charac- 
teristics (see fig. 4) could be input to the numerical solution for any given 
q_,aj combination, only small differences resulting in the shield motion. By 
calculating shield-half Inertias for a symmetrical configuration with the same 





Value of parameter for - 1 

Parameter 

Fire 

Long Nimbus 

XI, in. (m) 

1 ^ 5.20 

( 1 . 148 ) 

93.33 

( 2.371) 

X2j In. (m) ..... . 

50.10 

( 1 . 273 ) 

96.24 

( 2 . 444 ) 

in. (m) 

14.60 

( 0. 371) 

17.86 

( 0.454) 

y2, in. (m) 

13.90 

( 0 . 553 ) 

17.68 

( 0.449) 

leg lbf-sec 2 _ln. 

(N-sec^-m) 

560.00 

( 40.670) 

2964.00 

( 334.870) 

leg (D, Ibf-sec^-ln. 
(N-sec^-m) ...... 

504.00 

( 56.940) 

3324.00 

( 575.550) 

Weight ( 1 ), Ibf (N) . . 

156.00 

(695.920) 

304.77 

(1355.680) 

Weight 0 , Ibf (n) . . 

147.00 

(655.890) 

P86.09 

(1272.590) 

Ic, in. (m) 

24.50 

( 0.622) 

131.50 

( 3 . 340 ) 

It, in. (m) 

117.00 

( 2.972) 

224,00 

( 5.690) 

In, in. (m) 

10.05 

( 0.255) 

10.03 

( 0.255) 

Is, In. (m) 

76.00 

( 1 - 930 ) 

76.00 

( 1.930) 

Il3, in. (m) 

0.68 

( 0.017) 

0.68 

( 0.017) 

yg^, in. (m) 

12.25 

{ 0.311) 

12.25 

( 0 . 3 U) 

yag, in. (m) 

11.88 

( 0.302) 

11,88 

( 0.30a) 

bi, in. (m) ...... 

1.03 

( 0,026) 

1.03 

( 0 . 026 ) 

b 2 , in- (m) ...... 

1.00 

( 0.025) 

1.00 

( 0 , 025 ) 

d, in. (m) ...... 

62.43 

( 1.586) 

62.43 

( 1.586) 


Figure 4 . - Nominal shield physical data. 
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mass per -unit length characteristics as the short Nimhus or Fire shield^ deploy- 
ment trajectories were obtained for several combinations. After com- 

paring the motion of the near symmetrical halves with the motion of an assumed 
symmetrical half^ the agreement indicated that use of the input data for 
shield 1, the shield without the nose cap (fig. , was acceptable for rigorous 
motion analysis. A parallel investigation was also conducted for the long 
Nimbus shield without nose cap. 

Comparison with other work.- In lieu of experimental results to corroborate 
the analysis^ and in order to verify further the use of actual shield data in 
the analysis^ results for a O.lg vehicular acceleration case without aerodynamic 
loading for one Fire shield half were compared with the results of an inde- 
pendent analysis of reference k which uses a simultaneous solution to solve for 
the shield-separation dynamics. This comparison yielded the agreement indicated 
in figure 5. 

Aerodynamic Effects 

The prime purpose of this paper in including aerodynamic influence is to 
ascertain whether this parameter is important in the upper atmosphere and to 
determine whether there is sufficient justification for requiring more appropri- 
ate aerodynamic data for specific applications. None of the available published 
reports on shield deployment (such as refs, k and 5) treat the influence of 
aerodynamic forces. Although heat shields are generally deployed at very high 
altitudes^ it will be shown that significant influences on deployment trajecto- 
ries can be experienced even with very low dynamic pressures. In the applica- 
tion of this analysis where consideration is given to the importance of aerody- 
namic loading^ it should be observed that induced aerodynamic loading resulting 

from deployment veloc- 
ities (9^ y^ x) have 
not been considered; con- 
s e quen t ly y re suit s for 
cases in which aerody- 
namic loading is present 
are valid only when the 
shroud velocities are 
small compared with the 
velocity of the moving 
coordinate system. 

The aerodynamic 
characteristics of the 
complex time-vaiying 
geometry presented by the 
opening clamshell design 
are not amenable to ana- 
lytical solution. These 

Figure 5-- Comparison of analytical results for Fire or characteristics can only 

short Minhus shield without nose cap for a case of O.lg 
vehicular acceleration and zero dynamic pressure. 
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were obtained by extrapolation of wind-tunnel data (refs. 9 10 ) from models 

similar to the short Nimbus heat-shield halves used with Project Fire. These 
data do not include interference effects which might occur between shield 
halves during the early phases of opening and between a shield half and parent 
vehicle prior to movement away from the immediate vehicle region. Also, in 
establishing the aerodynamic characteristics, the aerodynamic center and the 
planform centroid were assiuned to be coincident for the shield half (ref. lO). 
In order to ensure the adequacy of the assumed aerodynamics used in the anal- 
ysis, individual variations in aerodynamic lift and drag, respectively, were 
examined for a representative q,a^ case. The results, as presented in the 
following table, indicate that the lift parameter variation of ±20 percent 
caused a corresponding clearance change of 6.9 percent and -7-7 percent. The 
drag parameter variation of ±20 percent caused a clearance change of only 
-I.98 percent and I.5I percent. Therefore, although changes in the lift param- 
eter have a greater effect on clearance, neither parameter variation affected 
the clearance distance appreciably. 


Aerodynamic varia- i 

tion, percent of 
nominal* 

Clearance 

L 

D 

Inches Meters 

138.046 3.506 

147.564 3.748 

127.412 3.236 

135.510 3.437 

i4o.134 3.559 

Nominal 

20 

-20 

Nominal 

Nominal 

Nominal 

Nominal 

Nominal 

20 

-20 


Percent clearance 
change from nominal 


6.9 

-7.7 

-1.98 

1.51 


*Nomlnal case: a2 = O.lgj q = I.65 lbf/ft2 (79,002 N/m2) 


Combinations 

Separation bounds . - In the following investigation, the aerodyTLamle q 
effect on shield deployment is considered with the aj loading since the two 
conditions occior in combination for actual launch vehicle trajectories. Fig- 
ure 8 Indicates the limiting q,a2 combinations under which separation can 
theoretically be effected for the Fire and long Nimbus shields by using the 
nominal spring actuator characteristics. The bo-undary line outside of which 
separation theoretically cannot occur is established by whether 0 vanishes 
during phase 2 or phase The spring locations are shown in figure 4, and 
the characteristics are given in the following table: 


Spring 

k 

Remarks 

Ibf/ln. 

N/m 


Forward spring actuator 

15.644 

2739.67 

Same values for Fire 

Aft kickoff spring 

102.564 

17961.62 

and long Nimbus 
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Results (fig. 8 ) show that at zero q, the long 
shield can be separated at a greater vehicle 
3.1 level than the Fire shield. Two primary 
reasons for this anomaly are the facts that (l) 
the long shroud center of gravity can be driven 
over "top dead center" (9 + 9co = 90^) "by 
rotating the shield through a smaller angle 
than is possible for the shorter shroud, and 
( 2 ) the inertial characteristics of the longer 
shield in rotation are more favorable than 
those for the shorter shield in sustaining 
motion once motion is initiated, (See fig. 4.) 
The larger aerodynamic influence on the long 
Nimbus shield accoimts for the increased sensi- 
tivity to dynamic pressure indicated by the 
greater negative slope of the curve for the 
long Nimbus shield in figure 8 . In estab- 
lishing separation bounds, notice also that in 
the initial phases of motion, the aerodynamic 
lift as well as drag oppose separation as shown 
by the negative lift force for 9 < 12,5^ 

(9 < 0.218 radian) in figure 6 . 

Atlas trajectory parameters ,- Since the 
Nimbus-type shield is commonly used on the 
Atlas launch vehicle, it is useful to present 
the possible flight times in a nominal Atlas 
trajectory during which the respective short 
and long shields can be deployed. These flight 
times can be ascertained by relatifig the lim- 
iting conditions from figure 8 to the actual 
trajectory conditions, the results being pre- 
sented in figure 9* Separation is theoreti- 
cally possible for a short time interval tp 



0 400 800 1200 ISoo 

Dynamic pressure, N/m^ 


Figure 8.- Locus of limiting 
conditions for separation. 
Fire and long Nimbus 
shields . 


immediately after launch. With the exception 

of possible applications to low-altitude abort systems, however, normal mission 
requirements will generally dictate shield separation subsequent to maximxam q. 
Although results indicate separation is feasible in the interval t2 fol- 
lowing BECO (Booster Engine Cutoff), it is obvious that the least complicated 
separation and greatest post separation clearance will be obtained in the non- 
thrusting region tj subsequent to SECO (Sustainer Engine Cutoff). 


Post Separation Clearances 

The clearance criterion .- The fact that either shield can be theoreti- 
cally deployed at any q^a^ combination within its respective envelope 

describing the limiting cases (fig. 8) does not necessarily imply that all 
these combinations are acceptable. For the possible separation conditions, 
another criterion to be satisfied is that the post separation trajectory of 
this type of shield satisfy a predetermined minimum safe clearance requirement. 
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Figiare 9.- Possible separation times during nominal Atlas trajectory. 

Long Nimbus and Fire shields. 

Shield trajectory.- Shield clearance trajectories are presented in fig- 
ures 10 and 11 for several assimed flight loading cases. For clarity, the 
shield tumbling motion has been shown for one case in each figure. The rota- 
tional velocity of the deployed shield yields a series of minima in the curve 
defining the clearance. These minima are defined as first minimum, second 
minimum, etc., as they progress downstream in relation to the parent vehicle. 
These typical points are indicated by the circular symbols on the 

q = 16.5 lbf/ft2 (790.025 N/m2) curve of figure 10. The aj values shown 
thereon indicate the vehicular inertial accelerations. 



Figure 10.- laical near point clearance trajectories. Fire shield 

with Atlas vehicle. 



itoo, 


200 100 


5-0 2.5 



-2.5 -5.0 -7.5 -10.0 

X dlQplaf'ement, m 


-12.5 -15.0 -17.5 -20.0 -22.5 


Figure 11.- 0?ypical near point clearance trajectories for 0.4g accel- 
eration with varying aerodynamic loads. Long Nimbus shield with 
Atlas vehicle. 


Comparison of clearance data for conditions.- The clearances as 

described graphically in figures 10 and 11 are presented in figures 12 and 13 y 
respectively, for a number of separation-possible conditions. The rea- 
son for the horizontally asymptotic behavior of the low q curves is that the 
separation angle for both long and short shields under low q, low a^ condi- 


tions causes shield post separation trajectories which are nearly perpendicular 
to the vehicle flight path. It is interesting to note from figure 12 that for 
all nonzero q cases shown, a constant first minimum clearance is indicated at 
the point of 0.45g vehicular acceleration. A similar condition exists for the 
long Nimbus shield at a vehicle acceleration of O.lg. Since the minimum 
clearance specified by the booster contractor for Project Fire was 120 inches 
(3«048 m), it is apparent from figure 12 that this clearance can only be 
obtained for very low q^a^ com- 


binations. From figure 13, it is 
seen that when the longer shield ^ 

is used, a greater latitude in t 

parameter selection will provide 
the prescribed 120-inch (3.048 m) ^ 

clearance. | 

g 

Constant-clearance loci . - " 

■ ■ — ■ a) 

Families of curves for common ^ 

first minimum clearances have been ^ 

extracted from data such as those 
given in figures 12 and 13 and are 
presented in figures l4 and I 5 for 
the short and long shields, 
respectively. The dotted curve of 
figure ih corresponds to the point 
of common intersection noted for 
the curve families of figure 12. 

The corresponding condition for 



0 .5 I'lo ' 1.5 " 2.0 ' 2.5 '* 5^0 

Clearance, m 

Figure 12.- First minimum clearance. Fire 
shield to Atlas vehicle . 
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1.5 2.0 2.5 5.0 3.5 '•-.O 4.5 5-0 5.5 

Clearonce , m 


Figure 13 •- First minimum clearance. Long Nimbus shield to 
Atlas vehicle. 



0 


4oo 800 

Dynonilc presGure, 


1200 

H/m® 


the longer shield is given by the 220-inch 
(5*588 m) clearance curve of figure 15- 
Notice^ however^ that figure l4 indicates that 
as q increases from zero to approximately 
3.5 lbf/ft2 (158^00 ir/m^) for a given g load^ 
post separation clearance is significantly 
decreased. For any a^ level greater than 
about 0«^5g^ a q- increase above 3*5 Ibf/ft^ 
(158.00 N/m2) becomes a favorable factor for 
increased clearance. Also^ for a^ < 0.^5g^ 

first minimum clearance tends to decrease with 
increase in dynamic pressure. To obtain maxi- 
mum clearance^ the results indicate that sep- 
aration for the short shield should be effected 
as near as possible to the zero q^ zero a^ 


condition. It should be noted that the 
120-inch clearance^ as required for the Project 
Fire mission^ can only be obtained with the 
highly restrictive combinations con- 

fined in the limited lower left-hand region of 
figure 14. Although the clearance family of c-urves for the long Nimbus shield 
is similar in shape to the family of curves for the short shield, it is seen 
from figures l4 and 15 that for the same q^a^ condition, the clearance values 


Figure 14.- Loci of common near 
point first minimimi clearance 
for separation-possible con- 
ditions • Fire shield. 


are conslderalily greater for the longer shield. Figure 15 also indicates that 
a 120-inch (3 .048 m) clearance reqixLrement can he satisfied with a wider varia- 
tion of q.,aj comhinations for the longer shield, ranging up to 10.2 Ibf/ft^ 
(488.4 N/ni2) and 0.68g. The analytical miniinTJin of the fixed clearance loci 
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apparent in figures 1 ^ and 15 repre- 
sents the least desirable combina- 

tion under which separation can be 
effected with respect to clearance. For 
the Fire shield^ this point occurs con- 
sistently for a dynamic pressure between 
3 and 4 l-bf/ft 2 (l 43.6 and I9I.5 N/m 2 ) . 

For the long Nimbus shield, the trace of 
analytical minima varies from zero to 
al 30 ut 2 lbf/ft 2 (95.8 N/m 2 ). It is Impor- 
tant here to note that in the region to 
the ri^t of the bounding curves for 
figures l 4 and 15, the shield would fail 
to eject. However, a discontinuoiis rela- 
tionship exists here in that just to the 
left of the bounding cuarves, separation 
will occTir and will yield the clearances 
indicated. In actual flights, of course, 
sTifficient margins on both q and a^ 
should be considered to allow for design 
parameter Tincertainties . In figure 15 ^ 
to the right of the trace of the minima, 
an increase in q causes a significant 
increase in clearance. It appears that 
aerodynamic loading in this q^a^ region 
is very favorable to post separation 
clearance . 


Spring Energy Effects 



0 100 200 300 koo 500 6 oo 

Dynamic pressure, K/m^ 

Figure 15 •- Loci of* common near point 
first minimum clearance for 
separation-possible condition. Long 
Nimbus shield. Some representative 
metric values of clearance are given 
by quantities in parenthesis. 


The Nimbus shield configuration con- 
tains spring systems which are normally 
considered to be fixed in size and posi- 
tion in their respective components. For example, the forward springs, or 
thrusters, are located at the same longitudinal position relative to the shield 
nose regardless of shield length. Similarly, the aft kickoff springs are always 
located at the shield separation ring (fig. i). Both forward and aft spring 
mechanisms are designed so that, for symmetrical halves, the spring forces 
always act in a direction perpendiciilar to the vehicle longitudinal axis. 


An investigation to determine the effect of variations in spring constants 
on shield deployment yielded the q^a^ results indicated in figures l6 and 17 
for the short and long shields, respectively, with the 100 percent or nominal 
spring curves replotted from figure 8. The curves describe the limiting condi- 
tions under which separation can be effected for the spring variations consid- 
ered. It should again be noted that not all the possible separation conditions 
to the left of a bounding curve imply acceptable separation conditions in view 
of minimum clearance requirements. 
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Upper variation of the spring constants was kept within reasonable limits, 
as dictated by the mechanical configuration of the existing Nimbus shield. For 
the spring variations considered, the same percent change was applied to for- 


ward and aft springs. It should be 
noted, as is evident from figures l 6 
and 17, that appreciable increases in 



0 56b 1000 1500 2000 2500 

Dynatalc pressure, , K/ra^ 


Figiire I6.- Spring variation effect on lim- 
iting separation condition. Fire shield. 
Separation not possible in region to 
right of boundary for a given spring 
variation . 



Dynamic pressure, q , N/m^ 


Figure IT*- Spring variation effect on 
limiting separation conditions. 

Long Nimbus shield. Separation not 
possible in region to right of 
bomdary for a given spring 
variation . 


the magnitudes of separation-possible 
q,a^ combinations are obtained through 

increasing the nominal spring charac- 
teristics. Results of the spring- 
variation effect on minimum clearance 
of the Fire shield to the Atlas booster 
vehicle for a representative 
8.25 Il3f/ft2 (395.01 N/m2), 0.5g con- 
dition are plotted in figure 18 . A 
similar curve is presented for the 
long Nimbus shield in figure I9 at the 
8.25 lt)f/ft2 (395.01 W/m2), 0.5g 
condition. 

Even though appreciable changes 
in separation bounds were noted in 
figures I6 and I7 with changes in 
spring rates, it is apparent from the 
comparative curves of figures I8 and 
19 that surprisingly little change in 
post separation trajectories would be 
obtained. 


Shield Length Effects 

In the foregoing discussion, no 
direct comparisons of shield length 
effect on deployment have been pre- 
sented. Since the Nimbus shield con- 
figuration is employed in many appli- 
cations and proposals with length 
variation of the aft cylindrical por- 
tion as the only significant geometiy 
change, two representative lengths 
have been selected and used throughout 
this investigation. The Fire shield 
is essentially a shortened version of 
the shield used on the Nimbus vehicle. 
The Nimbus vehicle shield is referred 
to as the long Nimbus shield. Some 
interesting results are obtained in 
comparing deployment for the two 
lengths with otherwise like structural 
and q.,a2 conditions. Figure 20 
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5.0 2.5 0 -2.5 -5.0 -7-5 -10.0 -12.5 -15-0 -17-5 -20. O -22.5 

X displacement, m 


Figure 18.- Effect of spring variation on Fire shield near point trajec- 
tory* for given loading condition. 



Figure 19*- Effect of spring variation on long Nimbus shield near point 
trajectory for a given q^a^ loading condition. 

illustrates the wide difference in pitch rate between long and short shields 
for a representative q^a^ loading condition. Moreover, the advantages of the 
slow pitch rate are apparent in that the first minimimi clearance is delayed 
until the longitudinal displacement is well beyond the parent vehicle. The 
unfavorable phasing of the pitch attitude of the short shield materially affects 
its minimum clearance (as is seen in fig. 20) at the downstream displacement of 
approximately 500 inches (7*^2 m). For the inputs used in this phase of the 
analysis, the comparison of figure 20 shows that the long-shield deployment 
provides greater clearance than the short shield even though the two internal 
ejection systems have identical potential energy. The favorable deployment 
characteristics of the long shield are attributed principally to the increased 
aerodynamic normal force that aids separation and to the favorable angular atti- 
tude of the shield as it traverses the length of the launch vehicle. The help- 
ful influence of the assumed aerodynamics was verified for the long Nimbus 
shield by comparing two cases with the same acceleration loading but with dif- 
ferent dynamic pressure as shown in figure 21. 
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Figure 20.- Comparison of shield- length effect on post separation near 

point trajectory. 
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X displacement, m 


Figure 21.- Comparison of aerodynamic effect on post separation trajec- 
tory for 0.5g acceleration. Long Nimbus shield. 

CONCLUDING REMAEKS 


An analytical investigation and parameter study have been made on a spring- 
actuated clamshell heat shield of the Nimbus type and the concluding remarks 
pertinent to the study are as follows; 

1. Results indicate that deployment of the Nimbus configuration heat 
shield, in the lengths studied, is feasible under significant aerodynamic- 
pressure— vehicular-acceleration (q,a^ ) combinations . 

2. Aerodynamic lift and drag loading can be of definite advantage in the 
low vehicle acceleration region for effecting a greater clearance in the post 
separation trajectory. 


30 



5. The aerodynamic effect is more pronounced for the longer shield because 
of the increased planform area or lift. 

4. The long Nimbus shield^ although more limited in the separation- 

possible loading envelope^ will attain a more desirable clearance profile 

than the shorter Fire shield for a given condition. 

5 . To obtain maximimi clearance when \ised on the Atlas launch vehicle in a 
nominal trajectory, both short and long Nimbus shields should be deployed after 
sustainer engine cutoff. 

6. Appreciable variation of the actuating spring rate affects the magni- 
tude of q,a 2 separation-possible conditions proportionally, but has negli- 
gible effect on the post separation trajectory. 

7 . The primary effect of shield length on deployment is in the shield rota- 
tional rate reduction at separation due to the increased rotational inertia 
associated with lengthening. 

The approximate aerodynamics used in this stuc3y appreciably influence 
shield separation conditions and trajectories when deploying in the detectable 
dynamic-pressure region. The significance of the aerodynamics indicates the 
importance of obtaining accurate wind-tiannel data on the aerodynamic charac- 
teristics of a clamshell heat-shield configuration for various stages of 
opening . 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 
Va., November 2, 


1965. 
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APPENDIX A 


NUMERICAL METHOD FOR SOLUTION OF THE PHASE 2 
ROTATION EQUATION OP MOTION 


The step-solution recurrence relations valid for solution of the initial 
value propagation problem, when linear acceleration in the interval is assmed 
and the time interval is considered as the independent variable, are as follows 


where 



(Al) 


(A2) 


(A3) 


and k represents the kth Iteration within a time step. Substitution of the 
expression for the third derivative into eqmtions (A1) and (A2), and sin 5 >lifi- 
cation yields 


. U) 


’n+1 


=n+l 


(k) 


^n+1 


(k) ^ 


0n+l^“' + 20. 


* 


At; + 0. 


n 


(a4) 

(A5) 


and from equation (24) 

l^ajA - rj^cos - cos(0n+i + 0^0 

+ k^|s - pjcos 0po - cos(9n+i + 9po^ 

“ [®n+l^^^ + ^2“]° cos 9co + Ln+i^^^c 




(a6) 
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APPENDIX A 


The method of solution is as follows ; By using the initial conditions of 
^ arbitrary value as an initial choice for ^n+l'" ' ^ equa- 
tion (A5) can be solved. Equation (a 6) is then solved by proper selection of 
L, and D as determined by next step is to obtain an 

Improved value of equation (a 6). The process is repeated until 

convergence in ©n+l obtained. By using this value, equation (a 4) is then 
solved. Also, at the end of each time increment, the following relations to 
deteimine whether the spring extensions are maximum must be evaluated: 


A - r|^cos 

0ro - cos (e^+i + 

0ro)] 

- ye 

^ ^1 

(AT) 

B - pj^cos 

0po - cos ( 0^+1 + 

9po)J 

- ye 

^ C2 

(A8) 


where and are prescribed positive limits ^==0. Since the arbitrary 

value of At selected for the problem may not permit these relations to be 
satisfied at given time increments to within the predetermined error boimds e , 
a method of subsequent halving of At is used as follows: If, in proceeding 

with the solution, a time increment is taken in which a change of sign occurs 
in the sum of the left-hand terms of relations (AT) or (A8), the solution 
has proceeded beyond the limit, and it is necessary to return to the time step 
prior to the occurrence of sign change, halve At, and proceed as before, with 
the new At. The sequence is continued until either of equations (A7) or (a 8) 
is satisfied, at which time the term involving that relation (ka or kb) van- 
ishes from equation (a 6). When both equations (A7) and (a 8) have been satis- 
fied in this manner, the step solution is then continued with the original At 
until translation on the pivot begins. 


55 



APPENDIX B 


NUMERICAL METHOD FOR SOLUTION OF THE PHASE 5 EQUATION OF MOTION 
DESCRIBING TRANSLATION ON PIVOT 


The noniterative recxirrence relations for the translational acceleration, 
displacement, and velocity on the pivot, with the assumption of linear accel- 
eration, are as follows; 



®n+lV ®n+l - ^*sin 


Ln+ll 


m 


(sin i|f - |ig cos \|f) 


rr. 


- \|0n+l(h 


cos 0n+i - x*sin ©n+l 




) - 


X cos 9n+i + tL sin 0. 


n+l> 


+ j(cos t + M'S 


(Bl) 


®n+l = (®n+l ®n)% + (®2) 

2 

®n+l (®n+l ®n 

Each time equation (a6) is solved for *0*n+l convergence, the res\ilting value 

is used to solve equations (a 4) and (A5) as before. Substitution of these 
quantities and the resTilting aerodynamic effects thus permit direct solution 
of equations (B1) to (B5) for each time interval. When the shroud has trans- 
lated on the pivot a distance e^, against which quantity equation (B5) is 
checked at the end of each time interval, shield deployment is effected. Notice 
that the interval halving scheme^ may again be req\iired in order to obtain 
em.2. = ®r within the prescribed error. 
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APPENDIX C 


NUMERICAL METHOD FOR SOLUTION OF THE PHASE h 
FREE-FLIGHT EQUATIONS OF MOTION 


By the s^e iterative method as employed in phase 2, convergence is now 
obtained on 0 for each time step by the following equations: 


^n+1 


(k) 



(Cl) 


%+l 


(k) 


(k) 


+ ©nAt + ©n 


ft ( k ) _ ^+1 

^n+1 “ 


(k) 


"eg 


(C2) 

(C5) 


from eq\iatlon (47) where k again denotes the kth Iteration within a time step. 
At convergence on for each At, selection of L and D at the corre- 

sponding angular displacement ©n+l permits direct solution of the following 
noniterative recurrence relations: 


from equation (^5) 


^n+1 = 


^+1 

m 


/•• , \At . • 

^n+1 “ (^n+1 ^nj 2 


and 


^n+l = (yn+1 + + y^At + y, 

Also 

W = - H(®n+1 + 


(C4) 


(C5) 


(c6) 


(C7) 
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APPENDIX C 

from eqijation (^6) , 

^+1 = (^+1 + (C8) 

and 

2 

%+l = (^+1 + 2Xn)^ + Xj^At + Xn (C9) 

The process is continued as long as required. The limit imposed could be some 
maximum value. of time or displacement as required. 
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”The aeronautical and space activities of the United States shall be 
conducted so as to co7itribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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